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Abstract 

The Darboux transformation applied recurrently on a Schroedinger operator generates what 
is called a dressing chain, or from a different point of view, a set of supersymmetric shape 
invariant potentials. The finite-gap potential theory is a special case of the chain. For the 
finite-gap case, the equations of the chain can be expressed as a time evolution of a Hamiltonian 
system. We apply Sklyanin's method of separation of variables to the chain. We show that the 
classical equation of the separation of variables is the Baxter T-Q relation after quantization. 
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1 Introduction 

What is the most universal method of solving completely integrable models? From Sklyanin's 
point of view || it is the separation of variables in its most general form. Therefore, it is desirable 
to understand old methods of integration in the light of the modern approach to the separation 

of variables. One of the most important "old" techniques applied with success to many physical 

* email: olipan@hsph.harvard.edu 

t email: crasinariu@popmail.colum.edu 



1 



and mathematical problems is the Darboux transformation. In short, the two main themes of 
this paper can be best described by two keywords: the Darboux transformation and the Sklyanin 
method of separation of variables. The Darboux transformation is, for example, a fundamental tool 
in the supersymmetric approach to quantum mechanics and in the theory of dressing chains GJ, 16]. 
Though the model is different from the DST model studied in |ij this paper and are closely 
related. The Sklyanin method aims to connect the separation of variables as we know it from the 
Hamiltonian mechanics with the new techniques of exactly solving mathematical physics problem, 
namely the Inverse Scattering Method and its quantum version ||. The worked example in this 
paper is the dressing chain representation of the finite-gap potential theory. We review the main 
ideas from the Darboux transformations with an emphasis on the Hamiltonian view-point on the 
finite-gap theory, following A. P. Veselov and A.B. Shabat GJ. We then work the Sklyanin method 



H> 111 on the finite-gap theory. We find the canonical separated variables for the dressing chain 
representation of the finite-gap potential theory. From the Quantum Inverse Scattering Theory 
point of view, the equation for the separated variables is the classical version of the Baxter T- 
Q equation. Then the quantum version of the finite-gap theory is presented, together with the 
corresponding R-matrix and the Baxter Q-operator. The conclusions and outlook will close the 
paper. 



2 The Darboux transformation and the Hamiltonian approach for 
the finite-gap theory 

Consider the Schroedinger operator for a potential u(x), TL = —D 2 + u(x) where D = d/dx. 
Factorize it product of two first order operators 

H = A* A , (2.1) 

where A = D — f(x) and A* = —D — f(x). The Darboux transformation sends 7i\ = A* A into 
TL2 = A A* + al where 1 is the identity operator and a is a constant. The functions f\ and 
obtained from the factorization 

H i = -{D + f i )(D-f i ) , (2.2) 
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are related by the equation 

(A + / 2 )' = /i 2 -/£ + «, (2-3) 
where the prime means differentiation with respect to x. From the supersymmetric quantum 



mechanics point of view, equation ( [2 .3; ) is exactly the shape invariance condition [16], written in 
terms of the superpotentials W% = —fi . 

We can continue this process of factorization and Darboux transformation and obtain a chain 
of equations 

(fi + fi+i)' = + * = 1,2,... . (2.4) 

This chain is called a dressing chain. In what follows we will consider only periodic chains 

fi = fi+N, OLi = Cti+N , (2.5) 

where the period N is a positive integer. Supersymmetrically speaking, periodic chains correspond 



to the cyclic shape invariant potentials [17]. The properties of the dressing chain depend drastically 
on the period N and the sum a = at\ + . . . + ajv- There are four cases to be considered depending 
if N is even or odd and a is equal to zero or not. The finite- gap theory, to be studied in this 
paper, corresponds to the case N odd and a = 0. In this case the chain is a completely integrable 
Hamiltonian system. We will use for the case a = instead of the constant Q!j, the constants 
related with a by — (3i+i = a%- If we regard the variable x in fi(x) as a time variable, then the 
dressing chain expresses the time evolution of the variables fi(x), i = 1 . . . N. Their evolution is 
generated by the Hamiltonian 

H = J2 N (If* + &fi) > ( 2 ' 6 ) 

with the Poisson bracket 

{fijj} = (-l)O--iWAT j (2?) 

{fiJi} = 0. (2.8) 

This Poisson bracket is not canonical. To obtain canonical variables, first let us denote by gi a set 
of variables defined by 

9i = fi + fi+i (2.9) 
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with the Poisson structure 

{^,^_i} = l. (2.10) 

all other brackets being zero. Though the variables gi seem to be redundant at this point, later on 
they will prove to be useful. Now, the canonical variables (JQ,Xj) 

{X t , X,} = { Xi , Xj } = , {Xi, Xj } = Sij , (2.11) 

will generate the Poisson structure for the variables g% if 

g i = X i + x i+1 . (2.12) 

As we mentioned before, the finite-gap case (N = 2n + l = odd) is completely integrable. Therefore 
the Hamiltonian responsible for the time evolution belongs to a set Hq,Hi, ...,H n of independent 
involutive Hamiltonians. To show this, we use the method of inverse scattering theory. The Lax 
matrix build on canonical variables (Xj,Xj) is 

L?\u)=( ^ 1 ) , (2.13) 

\ XiXi + (3i + u Xi J 

where (x) = (Xi,Xi) and u is a complex parameter. Then we construct the monodromy matrix 

L {x \u) = n Lf\u) = LP(u)L%Uu) . . . L^(u) (2-14) 

i=N 

and take its trace: 

TN (u) = TrL(u). (2.15) 
The trace tn{u) generates the set of involutive Hamiltonians 

t n (u) = H lU n + Hsu 11 - 1 + ... + H 2N+1 . (2.16) 

The Hamiltonian H3 is just the Hamiltonian for the chain. The fact that the set of the Hamiltonians 
is involutive, is a consequence of the classical r-matrix identity. Denoting by idi the unit 2x2 
matrix and introducing the notations for the tensor products = / <8> id,2, = id 2 <S> I, we have 

{lf\ Ul llf\u 2 )} = [n 2 (ui - U2 ),zf ) (« 1 )zf ) (u 2 )]^ , (2.17) 

where 

ri 2 («) = ~Pi2 , (2-18) 



and V12 is the permutation operator in C 2 <S> C 2 . It is important to notice that although the Lax 
matrix is written in terms of the variable (X, x), the Hamiltonians Hi ( |2.16| ) depend only on the 
variables gi, i = 1, 2, . . . , N. Moreover, it is possible to generate the Hamiltonians H{ as a trace of 
a monodromy matrix written directly in terms of the chain variables fi 

F(u) = l[l fi 1 ] , (2.19) 

i=i V fi+Pi + u fi J 

TN = TrF(u) . (2.20) 

The connection of the dressing chain with the finite-gap theory for the Schroedinger operators is 
described in [|]. For the finite-gap theory see also |fD| , |j| and the references therein. Here we 
will emphasize only those notions which will be important later on. To each solution fi(x) of the 
dressing chain, or in other words, to each solution of the time evolution of the Hamiltonian system, 
corresponds a sequence of N finite-gap Schroedinger operators 

Hi = -D 2 + Ui {x) , (2.21) 

where the potentials ui{x) are given by 

Ui = f'i+f? ■ (2.22) 
The spectral curve of these operators is the spectral curve of the monodromy matrix F{u) 

det(F(\) - fj, 1) = , (2.23) 
where 1 is the unit matrix, and it can be written as 

N 

/ u 2 -r 7V (A)/i-J](A + A) =0 . (2.24) 

i=i 

From the Darboux transformation we can obtain a recurrence relation between the logarithmic 
derivatives of the Bloch eigenfunctions. The Bloch eigenfunctions ipi of the operator Tii are given 

by 

Hi^i = (A + Pityi . (2.25) 
Due to the Darboux transformation, two successive Bloch eigenfunctions are connected through 

ij>i+i = {D- fi)^i . (2.26) 



This implies that for the logarithmic derivatives Xi = D In ifii we have the recurrence 

X t = fi + • (2.27) 

Ji + Xi+1 



3 Sklyanin method of separation of variables 



To understand the Sklyanin method |S|, gXj 11, gj, let us start with an old example. For an 
Hamiltonian of the form 

H = \{p\ + ..+pl) + -{u iq l + ... + u n ql) , (3.1) 

we notice that the variables are not only canonical, {pi,qj} = Sij but also separated i.e. each pair 
(pi(t),qi(t)) lies on the curve 

pf + ojiqf = const . (3-2) 

In general, let us consider a Hamiltonian system having d degrees of freedom and integrable in Li- 
ouville's sense. This means that it is given a 2d-dimensional symplectic manifold and d independent 
Hamiltonians Hi in involution 

{H i ,H j } = 0, i,j = l...d. (3.3) 
A system of canonical variables A = {Xi}f =1 and \i = {fn}f =i satisfying 

{Xi, Xj} = {/j,i, fXj} = , {Xi, pLj) = (3.4) 

will be called separated if there exists d relations of the form 

W j (X j , fJ , j ,H 1 ,...,H d ) = . (3.5) 

For the dressing chain the variables (X, x) are canonical but not separated. We can raise then the 
question of how to find canonical separated variables for the dressing chain. If the integrable system 
is solvable via the Inverse Scattering Method then we can use a method proposed by Sklyanin to 
find the transformation from the canonical variables (X, x) to the canonical separated variables 
(A,/i). The desired transformation will be obtained as a composition of Backlund transformations. 
The next section is devoted to the Backlund transformation for the dressing chain. 
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4 Backlund transformations 



Following Sklyanin, we need to find a canonical transformation (we will use also the name Backlund 
transformation) from the variables (X, x) to (Y,y). The important property is that the canonical 
transformation will depend on the spectral parameter A. This parameter A will allow us to find at 
the end the canonical separated variables. Being a canonical transformation, the Poisson structure 
and the set of Hamiltonians Hi must remain unchanged when expressed in the variable (Y, y). Since 
Lax matrix L(u) is a monodromy matrix 

L^{u) = L%\u)...Lf{u)Lf\u), (4.1) 

we can transform the Lax matrices at each site L * (u) 

Mi(u - X)Lf\u) = L^^Mi^u - A) , (4.2) 

because the trace tjv(m) is invariant due to M^{u—\)L^ x \u) = L^(u)M^{u—X). To keep the same 
Poisson structure we ask that the matrices Mi(u) should obey the same Poisson bracket ( 2.17| ) as 



Li{u) obeys. Practically, we first have to choose one out of many matrices which obeys the r-matrix 



Poisson bracket ( 2.17 ) and then be lucky enough to find that (|4.2| ) has a solution Y(X,x),y(X,x) 



for every spectral parameter u. The solution will depend on the parameter A which is exactly what 



we want. The technical way to solve (4.2) is quite interesting. The idea is to use, besides the phase 
spaces (X, x) and (Y, y), two more spaces: (S,s) and (T,t). The phase spaces (S, s) and (T, t) are 
auxiliary spaces, the Backlund transformation being between (X, x) and (Y, y). With the help of 
these auxiliary spaces, we can write a version of ([L^) as 

(« - X)L\ X \u) = L { y\u)M?\u - A) . (4.3) 

To go back to ( |4.2j) we simply need to impose the constraints 

U = , Ti = . (4.4) 

We apply now the above method to the dressing chain. The Lax matrix is 

Lf\u)=\ | • (1.5) 
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Here (x) stands for the pair of variables (X, x) and u is the spectral parameter. We choose the 
matrix M to be identical with L. Because in (14.31) the index i is the same on both sides, we can 



drop it and write the matrix equation as 

s 1 \ / x 1 

u-X + sS S ) \u + I3 + xX X 
The solution to the system is 




(4.6) 



S = -x + Z, (4.7) 

X = -s + ±±£, (4.8) 
t — x 

T = -y + ^, (4.9) 
t — x 

Y = + (4.10) 

We remark that £ is a free variable. This is a consequence of the fact that the conserved Hamilto- 
nians depend only on the combination X; L + Xi+%. The generating function F\(yt,xs) is 

F x (yt, xs) = y(t - £) - (x - £)« - (A + 0) ln(i - x) , (4.11) 

from which we get 

— f ■ 

To simplify the formulas we choose the £ = x. The generating function becomes 

F\ = y(t -x)-(\ + P) ln(i - x) . (4.14) 
The constrains tj = , T{ = S%-\ give: 

X, = - Si + A+ _ ft (4.15) 
which can be solved, in principle, for Sj. Then 

it + J/t+i = -Xi+i + Xi + Sj+i - Sj_i . (4.16) 
In terms of the old variables gi = X{ + Xj+i the transformation reads: 

^ = Zj , (4.17) 

Zi—l 



9i — 9i+l + z i-l — z i+l > 

where gi are the transformed variables and 



(4.18) 



(4.19) 



Note that ( 4.17 ) and ( 4.18 ) are just the canonical transformations we were looking for. To obtain 



the concrete form of these transformations we need to solve ( 4.17 ) for Zi, and then to use these 
values to find gi. Because for an arbitrary A equation (|4.17| ) cannot be explicitly solved, we leave 
the canonical transformation in an implicit form. However for special values of A, the canonical 
transformation can be explicitly solved, as we are going to exemplify in the section 10, eqs. (10.2), 
(10.3) and (10.4). 

The first equation in ( |4.17| ) is the discrete Riccati equation. It can be linearized with the help 
of the following change of variables: 

A+1 (4.20) 



ipi 



We obtain 



tp i+1 = ipigi + (A + 



i-i 



* = 0, ...,N . 



The periodic boundary condition zq = zn implies 



(4.21) 



(4.22) 



From ( 4.14 ), the generating function for the canonical transformation is 



Y; 



dxi 
dyi 



N 



N 



®\(y,x) = y y F x (yiSi-i\xiSi) = ^^y^s^i - - (A + pi) ln(sj_i 



(4.23) 
(4.24) 

(4.25) 



i=i 



Here we denote by x = (x±, . . . , xn). Finally in terms of Zi the generating function can be written 



MS 



N 



®\(y, x)=^2 yi(~ z i^) - ( A + A) ln (-^-i) 



(4.26) 
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5 Canonical transformations and Darboux factorization 



To decipher the meaning of the variables z% present in the canonical transformation ( 4,17 ) we will 
use the knowledge obtained from the Darboux method of factorization. First, from ( f4.17|) , find Zi-\ 

- Zi-i = • (5.1) 

9i - Zi 

Then, compare this result with formula ( 2.27j ) that gives the recurrence relation between two 
logarithmic derivatives of the Bloch eigenfunctions 



Xi = fi + 



fi + Xi+i 



(5.2) 



We obtain thus: 



(5.3) 



Zi-l = fi-Xi , 
which in terms of the superpotentials Wi = — fi reads as 

-Zi = Wi + Xi ■ 

Therefore the variable Zj taken with a minus sign, is the sum between the superpotential W» and 
the logarithmic derivative of the Bloch eigenfunction. It is interesting to obtain the time evolution 
of the variables Z{. From 



we get 



(A + I3i)ipi , 

f'i ~ 



H-l 



xl - fi + A + A 



(5.4) 
(5.5) 
(5.6) 

(5.7) 



6 Separated canonical variables 

At this point we have the canonical Backhand transformations (4.17). Let us use the symbol B\ for 



this transformation. Our goal is to find separated canonical variables. Following Sklyanin, we con- 
sider the composition B\ 1 ...\ N = B\ l o ... o B\ N of Backlund transformations and the corresponding 
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generating function F\ 1 ^_x N (y,x). If we treat A's as dynamical variables and y's as parameters then 
F\\...\ N {y^ x ) becomes the generating function of the ^-parametric canonical transformation from 
(X, x) to (fi, A) given by 

This transformation is not only canonical but also separates the variables. See || for details. 
Each pair (Aj, /ij) lie on a curve given implicitly by 

W(Xi,fXi)=0. (6.2) 

To find the curve W we use the parameter \i which is the variable conjugated to A 

= dX ' (6 ' 3) 

and search for a function /(/i) such that 

det(f(fi) - L(A)) = . (6.4) 

Then the spectral curve W is 

W(\i, f(fii); Hi) = det(f(^) - = . (6.5) 

For the dressing chain 

/(/i) = -e" . (6.6) 

To prove this, we will show that — e M is an eigenvalue for L(X) so the property det(f(fi) — L(X)) = 
is immediate. 



From the definition of fi ( |6.3| ) and from the generating function ( [4.26 ) we get 



z x ■ ■ ■ z n = -e" . (6.7) 

Now, by a simple computation 

L?>(A)| I =: ; _ J j ! . (6.8) 



This proves that for L = Ln ■ ■ ■ L\ the eigenvalue is zi • • • z^. 
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The spectral curve W(X,fi) can be expressed in terms of the trace tn and 0i 



N 

det{v - L(A)) = v 2 - t n (X) v + + bi) (6.9) 

i=l 

SO 

N 

e^ + r N (X)e^-H(X + b i ) = . (6.10) 

i=i 

Here 

t n (X) = H\X n + H^- 1 + ... + H 2N+1 , (6.11) 
where N = 2n + 1 and Hi, H 5 , . . . , H 2 n+i are integrals of the chain. 

7 Quantum case 

To get the quantum version of the theory described so far we will use the R-matrix approach. This 
will ensure the commutativity of the Hamiltonians Hi after quantization. Prom classical variables 
(x,X) we move to the quantum variables (x,d x ). The local quantum Lax matrix 

( x 1 \ 

L(u\x,d x )=\ (7.1) 

\ u + xd x d x J 

verifies the quantum commutation relation 

R X2 {ui - u 2 )L^( Ul )L^(u 2 ) = L^(u 2 )L^(ui)Ri 2 (ui - u 2 ) , (7.2) 

where 

Ri 2 (u) = u + Vi 2 (7.3) 

is the 5 , L(2)-invariant solution to the quantum Yang-Baxter equation [13]. The monodromy op- 
erator and its trace are defined like in the classical case. The commutativity of the Hamiltonians 
Hi 

[Hi,Hj]=0, (7.4) 

is a consequence of ( |7.2| ). The whole machinery of the Quantum Inverse Scattering Method can be 
put to work at this stage. We will limit to study only the Baxter Q-operator and the Baxter T-Q 
relation. The Q-operator will depend upon the spectral parameter A. Let us denote it by Q(X). 
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The interesting aspect is that the classical Backlund transformation B\ is the classical limit of the 
similarity transformation 

O^Q(X)OQ-\X) . (7.5) 
For details see ||. In the next section we will explicitly construct Q(X) as an integral operator. 

8 Q-operator 

For the Baxter Q-operator we require the three usual properties. First, it has to commute with the 
trace of the monodromy matrix t^{u) = nl=7v L(u\xi,d Xi ) 

Wu),Q(A)]=0, (8.1) 

second, it has to commute with itself 

[Q(X 1 ),Q(X 2 )]=0 , (8.2) 

and the last important property imposed is the Baxter T-Q equation, i.e. the Q-operator should 
satisfy a finite difference equation 

tjv(A)Q(A) = A(X)Q(X - 1) + B(X)Q(X + 1) (8.3) 

where A(X) and B(X) are two functions (not operators) of the spectral parameter A . 

We will follow |Q] and construct Q(X) as an integral operator 

N 

(Q(A)V)(x) = / dt dffYlRx+fii-iiUiXilti-iViMif) ■ (8-4) 
^ i=l 

Here dt = dtjsf-.dti and similar for alike symbols. If we introduce the 7^,-operator as 



(n x i>)(s,x) =JdyJ dtR x (s,x\t,y)^(y) , (8.5) 
the formula (|8.4| ) can be understand in the general sense of the trace of a monodromy matrix 

Q(A) = Tr tN n\ + p N _ x ...n^ Pl _ x . (8.6) 

In the notation R\(t, y\s, x) we recognize the auxiliary indexes s,t and the quantum indexes 
x, y. The Q-operator can be expressed as an integral operator 



(Q(A)V)(x) = / d yi ... I dy N Q x (x\y)^y) (8-7) 
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with the kernel 

Qx{x\y) = I dt N ... I dh Y[ R\+/3i-l(ti-i,Xi\ti,yi) ■ (8.8) 

J i=N 

After this general introduction, we move forward to find the concrete form of the operator 1Z\. The 
first property of the Baxter Q-operator, namely the commutation [tn(u),Q(\)] = is fulfilled if 
TZ\ is a solution of an equation similar to (177: 



M(u - X\s, d s )L(u\x, d x )TZ x = K x L(u\y, d y )M(u - X\t, d t ) , (8.9) 



where L(u\x, d x ) is the local quantum Lax matrix (7.1) and M(u — A) is another matrix which obeys 



the quantum commutation ([7.2j.) The main difficulty is how to chose the matrix M{u — A) so that, 
the equation ( |8.9| ) for 1Z{\) has a solution for every complex parameter u and by the other hand 
the Q-operator thus obtained has the required properties. The second property of the Q-operator 
comes from the Yang-Baxter equation which can be obtained from (|8.9| ) by a standard technique, 
see [13]. Returning to equation ( |8.9[ ) we take M to be of the same form as the Lax matrix ( [7.11 ). 



We obtain 




Rx(t,y\s,x) 



x 1 

u + xd x d x 

y i \ I i i 

u + ydy d y ) \ u- X + tdt d t 



R x (t,y\s,x)= (8.10) 



On the right hand side of the above equation, move R\(ty,sy) from the left side of the matrices 
product, to the right side. We have to change d x — » —d x and xd x — > — 1 — xd x . Then the equation 
becomes 

' s 1 \ / x 1 

u — A + sd s d s I \ u + xd x d x 

y i W t l 

u-l-ydy -d y ) \ u - A - 1 - td t -d t 
The solution is: 

R x (t,y\s,x) = Px 5(s - y)e^- x \t - x)'^ 1 . (8.12) 

We notice that R ~ exp(F\) for /? = 1. Due to the Dirac function, the solution is gauge 



R x (t,y\s,x)= (8.11) 



R\(t,y\s,x) 



independent, i.e. the solution does not depend on the free variable £ from ( [4.11|) , The 7^-operator 
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1|) becomes, after integration over y and changing the variable t — x = £ 

(K4,)(s, x)= Px j d^ C A ~V(z + £, s) , (8-13) 



or 



(Ki/j)(s, x) = Px s x / d^ r V(s + s" % s) . (8.14) 



The branch for the many valued function s A from ( 8.14j ) is fixed by making a cut along (— oo,0) 
and taking arg(s) G [— mir\. We have to specify the factor p\ and the integration contour (in the 
complex £ pane) in ( 8.1^ ). The integration contour is the Hankel contour for the Gamma function 



The previous formula inspired us to choose 



(0+) * 2m 

e^-*dg = . ( 8 .i5) 



PX = 2^ r(A + 1) • (8 - 16) 



Then 



(KxWKs, x) = ^-T{\ + l)s x f d^r A ~V(z + s-\, s) . (8.17) 

We are ready now to write the kernel of the Q-operator ( |8.8| ). From ( |3.12| ) and ( |8.6| ) we obtain 

N 

Q\{x\y) =Y[wi(\;y i - 1 ,y i ,Xi) (8.18) 

i=i 

where 

Wi(X;yi-i,yi,Xi) = -Lr(A + Aje*'*^-!-**)^ _ ^ )-a-A . (8 . 19) 
2m 

Therefore we have found an explicit form for the Baxter Q-operator ( |8,18| ). Next we are going 
to investigate the third property of the Q-operator, namely the Baxter T-Q equation. 



9 Baxter T-Q equation 



This last paragraph aims to show that the Baxter T-Q equation is the quantum version of the 
classical separation of variables ( |6.l[) . The computation parallels the one in Q. Start from the left 
side of the Baxter T-Q equation AO) 



[r(A)Q(AM(s) = TV 



/ dtdyl Y[ L(.M x i: d Xi) R X+l3i-l(ti,Xi\ti-l; 
J \i=N 



yi) Hv) 



(9.1) 
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We can integrate over ti and get 



[r(X)Q(X)i;](x)=Tr 



\i=N 



(9.2) 



where Wi are given by (8.19). 
Move all Wi to the left using 



L(\\xi, d Xi )wi = WiL(X\xi, d Xi ) 



with 



L(X\xi,d Xi ) 



or 



L(X\xi,d Xi ) 
At this point we can write 



Xi 1 

A + fa + Xid Xi In Wi d Xi In Wi 

Xi 1 

(X+Pihi-! , (A+ft) 



(9.3) 



(9.4) 



(9.5) 



[r(A)Q(AM(aO = / [] t^Tr [L(X\x N , d XN ) . . . L(X\ Xl , d xi )j ^{y) . (9.6) 

J i=N 

The last step is to perform a gauge transformation which leaves the trace invariant and make the 
matrices L(X\x{, d Xi ) triangular, so the trace will be easy to compute 



(9.7) 



With the help of the following gauge matrix 



N; 



1 



(9.8) 



the triangular form for L(X\xi,d Xi ) is 



Nr'-LiXlxud^Ni-i 



■(yi-i-Xi) 1 







(9.9) 



Hi-l—Xi 

The entries of the previous matrix can be expressed in terms of the Wi 

Wi(X-l) 



(yi-i - Xi) = -(A + A)" 
X + pi Wi(X + l) 



Wi(X) 



yi-i-Xi Wi(X) 
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(9.10) 
(9.11) 



so we get for the trace 



Tr (L(\\x N , d XN ) . . . L(X\ Xl ,d xl )) = JJ -(A + + II ^^^7^ ' ^ 

i=l 1 i=l 

The last result implies the Baxter T-Q equation 

AT 

r(A)Q(A) = - H(\ + Pi)Q(\ - 1) + Q(A + 1) . (9.13) 



Compare ( |9 . 1 3| ) with the classical result (6.10) written in the form 



N 



e» + r N - + bi)e^ = . (9.14) 



The connection is obvious if we quantify the canonical pair (//, A) as 



M -> ^ , A - A . (9.15) 



Then ( p. 14 ) becomes an operator acting on the Q-operator 

{e» + t n - + h)e-^j Q(A) = , (9.16) 
which is the T-Q equation up to a minus sign. To obtain a T-Q equation which exactly matches 



the classical formula, we have to chose for p\ the one in ( 8.16 ) multiplied with (— 1) . 



10 The case N=3 

It is instructive to study the case N = 3 which corresponds to the one-gap potentials. The trace 
( 2.15| ) of the monodromy matrix ( |2.14| ) for N = 3 is 



T 3( u ) = (ffl + 92 + 9z)u + g\Qi9z + 9103 + 92 01 + 9302 ■ (10.1) 

The variables §i, i = 1,2,3 are (|2.12| ): g\ = X\ + %2 , 92 = X2 + £3 and §3 = X3 + x±. The 
transformation -B(A) ( [4.17 , 4.1S| ) can be explicitly found for A = —0i, i = 1,2,3. For example, for 



A = — 02 we get 



.03 - 02 , ln „x 

9i = 93 H , (10.2) 

92 

03 - 02 , 

52 = ft " — + ^H' (10 ' 3) 

A -02 , ln ,s 

53 = 52 " ' (10 ' 4) 
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This transformation can be recovered from the Backlund transformations Tk, k = 1,2,3 from 
@, 0. Recall that T k is given by 

Tk(9k±i) = 9k±i ± — — , (10.5) 

9k 

T k (f3 k ) = Pk+i , (10.6) 

W fc+ i) = ft , (10.7) 

the remaining ft and gj being not changed. We also need to introduce the shift S acting as 

S{f3i) = ft_i , (10.8) 
S{ gi ) = 9i-i • (10.9) 

In terms of these last transformations, we can write 

B(-(3 2 ) = T 2 ST 1 . (10.10) 

We cannot recover T k from B(X) because of the difference in nature between these transformations. 



T k transforms the parameters ft so it changes solutions of one system of equations (2^) to solutions 
of another system of the same type ( |2 . 4j ) . The transformations B(X) change the solutions of the 
same system ( |2.4| ) among themselves. In this respect B(X) is an auto-Backlund transformation. 

We can try to solve the discrete Riccati equation ( [4.20 ) for ipi. In this case, we will get ipi,---,ip3 
in terms of ipo and ^4 

^i = | [V>4 - M9293 + A + ft)(A + ft)] , (10.11) 
^2 = - [V>4Si + V>o93(A + ft)(A + ft)] , (10.12) 
^3 = \ [M9X92 + A + ft) - V>o(A + ft)(A + ft)(A + ft)] , (10.13) 

with Z = 9x9293 + 5i(A + ft) + ^(A + ft). It obvious that Zj = ipi + i/ipi will depend on ip^ and V'o 
only through their ratio ipt/ipo. This ratio is not a free parameter because the periodic boundary 
condition zq = Z3 imposes a restriction on it. 

Though the r -functions are not one of the major players of this paper, it is worthwhile to 
mention the connection it has with canonical transformations. 
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The r -functions for the dressing chain were reported in [18] for N = 3. For N = 3 case, the 



dressing chain (2.4) written in variables gi, ( |2.9| ), is: 

-9o(9i ~ 92) + A) - 01 > 
-51(52 -0o) +P1-P2 , 
-92(90 -91) + P2- Po • 



9'x 



(10.14) 



92 



With the change of variables 



90 


= F[ 


-F2 + C, 




91 


= H 


- ^0 + c , 


(10.15) 


92 


= H 


- F{ + c , 





where the constant c is given by 3 c = go + g\ + 5% the system of equations ( |l0.14j ) transforms into 

F^ + F{' + (F^-F[) 2 -c(F^-F{)+p 1 = 0, 

F\ + F'-l + (F{ - F2) 2 - c(F{ - Ffi + fa = , (10.16) 
Fi' + FZ + W-Ftf-cM-Fti + Po = 0. 

The r -functions To, ri, Ti are now given by 

F = log r , Fi = logn , F 2 = logT 2 . (10.17) 

In terms of the r -functions, the dressing chain becomes a Hirota type system of equations: 

(L%-cD x + Pi)Tb-n = , 

(D 2 X - cD x + (5 2 ) n • t 2 = , (10.18) 
(Dl - cD x + /3 ) t 2 • t = , 

where for a polynomial P, the operator P(D X ) is defined as 

P(Ly F(x) • G(x) = P(d y ) F(x + y)G(x - y)\ y=0 . (10.19) 



At this point the goal is to work the canonical transformation ( 4.17 ) and ( |4.18| ) in terms of the 
t -functions. The canonical transformation ( 4.17[ ), ( |4.18] ) namely 



9i = Zi 



Zi-l 



(10.20) 
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and 



9% — 9i+l + z i-i ~ z i+l : 



(10.21) 



cannot be written explicitly as a formula which comprise only g~i and gi. Therefore we use the 
following strategy: given gi, we have to solve for z% in ( 10.20[ ), and then obtain the transformed 



variables cji with the aid of ( 10.21 ). As a result, the canonical transformation for the r -functions 
will be written in terms of the variables Zi. 



First, by simple manipulations of (10.14) and ( |10.15| ) we get 

p// 2 00+01- 02 

F = gm - c 



and similarly for F" and F%. Using (10.17) we get 

(log e -70 *' 2 t ) = 5-15-2 , 



(10.22) 



(10.23) 



1 ( c 2 + §p±fh-fh 



and similarly for t\ and T2- Here 70 = 

After we apply the canonical transformation we obtain the function tq given by 

(log e" 70 *' 2 fb) = M2 ■ 



(10.24) 



Now it is easy to write the canonical transformation for tq,t±, ti in terms of the variable 2, . Use 

X + 02 



(110.201) and (110.21Q in (|10.23|) and (|10.24|) and get: 

X + 0i 



loge'^ T 



loge 71 ti 



(loge" 72 ^)' 



zi 



Z2 



and 



(loge" 71 ^)" 



Z\ 

Z-2 



zo 

X + 02 
Z\ 

X + 0o 

Z2 



X + 02 
Z\ 

A + /3 

Z2 

X + 01 

zo 



Z2 

zo 

Z\ 



Zl 
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zo 



with 



I f 2 , 0O + 01-02 \ 1 ( 2 , 01 + 02 ~ 00 

70 = -2l C + 3 ' 71 = "2 C + 3 



Zl 

X + O 

Z2 
X + 01 

zo 



X + 0o 

Z2 

X + 0i 

zo 

X + 02 

Zl 



72 



(10.25) 



(10.26) 



1 / 2 , 02 + 00 - 01 



ff + 
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The meaning of the transformations as presented in ( |10.25| ) and (10.26) is that first we must 
factorize like in ( 10.25| ) and then use the factorization variables z% to obtain the transformed r 
-functions, like in ( 10,26| ). 



For the case N = 3, the canonical transformation must be carried on three times, each time 
with another A in order to obtain separated canonical variables, see ( |6.1[ ). In variables gi, these 
transformations read as 

Ai + A 



9i = 9i+i + Zi-i - z i+ i , gi = Zi 



9i = 9i+i + ~ Zi+1 , 9i = Zi 



9i + \ + Zi-i — Zi+i , 9i — Zi 



Zi-l 

A 2 + A 

Zi-l 

h + Pi 

Zi-l 



(10.27) 



If we add fj,%, fi2, /^3 given by (6.7), 



,A<1 



-ZlZ 2 Z 3 



-Z!Z 2 Z 3 



-ZlZ 2 Z 3 



(10.28) 



the variables (Ai, A2, A3, /ii, ^2, /X3) are separated canonical variables. The time evolution of these 
variables, inherited from the dressing chain ( |10.14| ) is such that, ( 6.10| ), 

e 2m + e w _ (Al + /3l)(Al + / 3 2 )( Al + /3 3 ) = , (10.29) 
and similarly for the pairs (A2, H%) and {X^,^)- Here Tjv(Ai) is given by (10.1) with Ai instead of 



u. 



If we use the notation = r , = f , = f and similarly for Zi, the transformations 



(10.27) for the r -functions, can be compactly expressed as follows: 
factorization: 

~t,x\{k)\" _ (m A fc + A+i ) (m A fc + p i+2 



loge"^ t 



"t+l 



"i+2 



(fc) 
Z i+l 



transformation: 



(lo ge --Vf«»)"=( z f'-^±|i 



i+2 \ ( (fc) _ Afc + $ 
(fe) 

i+2 



Z 



Here i = 0, 1, 2 modulo 3, and k = 1,2, 3. The variables are given by (gaJD as before. 
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11 Hamiltonian flow for the case N = 3 



Here we discuss the time dependence of the canonical variables (Xi, Xi). We regard the variable x 
(i.e. the space variable) in the system of equations ( |2.4D as being a time t for the Hamiltonian flow. 
The Hamiltonian that governs the motion in time is in variables /j ( p.9[ ) given by 

h = ^(A 3 + /! + /| + Ihfi + fah + Ms) ■ (11.1) 

It is useful to list all the variables which appeared so far 

91 = fl + h > fl = 2(51 - 52 + 93) , 51 = -Xi + %2 , (H-2) 

92 = f2 + h , f2 = ^(52 - 53 + 5l) , 52 = ^2 + X 3 , (11.3) 
53 = /3 + fl , h = ^(53 - 51 + 52) , 53 = ^3 + Xi . (11.4) 

The evolution of the variable x\ in time is given by 

dx 



J- = { X1 ,H}. (11.5) 



The Poisson bracket can be computed for two arbitrary functions / and g from 

df dg 



a,/? 

where (yx, ■ ■ ■ , 52iv) = 0-^ij • • • > ^iv> x i> • • • > ^tv)- I n this way we arrive at 

2^ = (A 2 + fl - f! + P1+P2- P 3 ){xi,X 1 } , (11.7) 

which gives the evolution for x\ knowing that the variables (X\,x\) are canonical, i.e. 

{x 1 ,X 1 } = l. (11.8) 

From the paper of Veselov and Shabat Q we even know the solutions for the dressing chain in 
terms of the elliptic Weierstrass "P-functions: 

2 V(t + di) - V{bi) 

We can integrate ( |11.7|) to 

2 Xl {t) = £(t + a 3 + 63) + £(* + a 3 ) + £(6 3 ) - Z(t + 01 + 61) - f (t + ai) - - 

e(t + a 2 + 6 2 ) - c ( t + 02) - e(^) + a + ih - & • (11.10) 
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Here £' = —V, bi = V{fii) and Oj+i — aj = 6j . 

In the quantum case the Hamiltonians are generated from the Lax matrix ( |7.1[ ). The trace of 
the monodromy matrix is: 

r 3 (u) = H lU + H 3 (11.11) 

where 

iTi = xi + x 2 + x 3 + d Xl + 9 X2 + <9 X3 , (11.12) 
Hz = d Xl d X2 d X3 +x 1 d Xl d X3 + x 2 d X2 d Xl + x 3 d X3 d X2 

+ x 2 x 1 d Xl + x 3 x 2 d X2 + xix 3 d X3 + x 1 x 2 x 3 . (11.13) 

An interesting feature of the quantum case is the absence of evolution of the variables = 
d Xi + Xj + i due to 

[ 9i ,H 3 } = 0. (11.14) 

This means that there is no quantum analog of the dressing chain in terms of the variables gi. 
The only surviving variables are Xj = d Xi and Xi which each separately evolve in time under the 
Hamiltonian ^3. 



12 Conclusions 

For the finite-gap potentials we have shown that the Darboux transformations can be viewed as 
canonical transformations if we apply the method of separation of variables proposed by Sklyanin. 
Not only the finite-gap case (N=odd, a = ) is interesting but also the other cases. For example 
the spectrum of the potentials which are solution of the chain in the case ct / and arbitrary 
N is described as following. The ground state is at zero energy; the next (p — 1) eigenvalues are 
Ei = Ylk=o a k > ' = 0) 1' • • • j (p ~ 2) 1 and all other eigenvalues are obtained by adding arbitrary 
multiples of the quantity a = «o + a i + • • • + «p-i- The general formula for the excited energy 
levels is: [17] 



1 

na + Y,ak ; {n = 0, 1, 2, . . . , 00 ; I = 0, 1, . . . , (p - 1)} . (12.1) 

fc=o 

The above potentials (also called cyclic shape invariant potentials) are a direct generalization of 
the harmonic oscillator. For N = 3 the potentials are Painleve transcendents 
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In this work we have also shown that there exists a quantum version of the dressing chain, 
namely the time evolution of the variables (X, x) under the Hamiltonian H\ . From here there are 
many ways to proceed. One way is along Bethe-Ansatz procedure. It will be interesting to find 
the spectrum of the quantum Hamiltonians. Also, each Hamiltonian Hi has its own time ti for 
evolution, so there must be a r-function which depends on all time variables r(ti,t2, ...,t^). The 
r-function for the dressing chain was reported in [^] for N = 3 in connection with the Painleve 



equations. See also [19]. 



We have analyzed the role of r -functions in connection with the factorization method and 
canonical transformation. The r -functions are an important tool in understanding integrable sys- 
tems and we believe that the connections between the r -functions and the Darboux transformations 
is worth to be studied further. 

Finally, a word about KdV. The finite-gap potential theory provides solutions of the periodic 
boundary problems for KdV equation. The KdV equation 

tlig^-fc — ^-^XXX ^^-^^-^XX ^12.2^ 

is a partial differential equation in two variables. One variable x we interpreted as a time variable 
associated with the Hamiltonian H\. To what Hamiltonian is the second variable, i.e. t, associated? 
For Pi = the Hamiltonian is H§ — H\H^. This result is buried in the paper [Q. What is the 
meaning of the Sklyanin separation of variables for KdV equation, both classical and quantum? 
From the conformal field point of view, the quantum KdV and the T-Q relation was already studied 
in the paper of [p^]. 

In conclusion, the Darboux transformation together with the new approach of the separation 
of variables is a promising research direction. 
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